On strategies for determination and characterization of the underlying event 
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^s^j ■ We discuss the problem of the separation and description of the underlying event (UE) within 

two existing approaches to UE measurement: the "traditional" method, widely used at Teva- 
tron, and a recently proposed jet-area/median method. A simple toy model of UE is developed 
i-C ' in order to understand how these approaches perform. We find that both methods are com- 

parably good for measuring average properties of the UE but the jet-area/median approach is 
O ■ favorable for determining fluctuations. We also use the latter method to study the UE from 

Q_) ' several existing Monte Carlo generator tunes. We investigate which characteristics of the un- 

, derlying event might be useful to measure in order to improve understanding of its properties 

and to simulate it well. These include transverse momentum density per area, intra- and 
inter-event fluctuations and correlations. 
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The hard processes studied at hadron collider are nearly always accompanied by the underlying 
event (UE). This predominantly soft activity affects a wide variety of high-pj measurements, 
e.g. by introducing a bias or by degrading kinematic jet reconstruction. Therefore, a good 
understanding and precise theoretical control over the underlying event is of great importance 
in order to fully exploit the potential of LHC. 
^ Study of the UE, both experimental and theoretical, faces, however, a number of problems. 

The first of them appears already at the level of definition since the very concept of the un- 
derlying event is ambiguous. This is because none of the events measured at hadron colliders 
has a clear boundary between the hard part and the UE. Also modeling of the UE is difficult. 
Most successful and widely used Monte Carlo models generate the majority of the underlying 
activity via the mechanism of multiple parton interactions. There are however other conceivable 
mechanisms which could contribute to the UE and which have not been given as thorough study. 
Questions include e.g. the role of correlations as well as the effects of possible contributions from 
the BFKL type radiation. 

Given the above difficulties, one can ask the question if, at least, it would be possible to 
measure the UE in, admittedly always to some extent arbitrary, but well defined and well 
controlled way. 

This leads us to addressing the following two questions:^ What do we really measure with 
existing methods of UE determination and which observables are interesting to measure? To 
answer the first question we develop a simple toy model and use it as a testing ground for two 
existing methods of UE measurement. To address the second question we study the UE from 
Monte Carlo models and identify a set of quantities which could serve as valuable characteristics 
of the UE. 



2 Relevant characteristics and measuring methods of the UE 



We concentrate on quantities related to energy flow. Those involve the main observable called p, 
which is defined as the amount of transverse momentum of UE per unit area. We will be also 
interested in rapidity dependence of p, its point-to-point fluctuations within a single event, a, 
and its fluctuations from event to event as well as the point-to-point correlations. 

Tw o m ethods exist which allow one to study the underlying event: the traditional ap- 
proach, 1^121 widely used at Tevatron, and the more recent area/median based approach. 
Certain freedom is present in each of the methods. Below, we describe the version used in our 
study. 

In the traditional method all particles within a certain rapidity range are first clustered into 
jets. Then, four regions in the transverse plane are defined based on the position of the leading 
jet, whose direction defines = 0: the "towards" region, < ir/3, an away region, 2-7r/3 < \4>\ < 
7r, and two transverse regions. The characteristic pt of the UE is defined as pt in one or both of 
those transverse regions. To further reduce the contamination from the perturbative radiation, 
the transverse regions are labeled, on an event-by-event basis, as TransMin and TransMax 
depending on their relative value of pt- The contribution of perturbative radiation to the average 
Pt in TransMin should be suppressed by an additional a s with respect to the average result from 
both regions TransAv. 

The area/median method is jet-b ased, exploits the concept of jet areas^ and can be carried 
out using the Fast Jet package. It starts by adding a dense set of infinitely soft partic les, 
ghosts, to an event. Then, all particles (real and ghosts) are clustered with the C/A algorithup^SI 
leading to a set of jets ranging from hard to soft. The typical UE scale in the event is defined as 



where Aj is a jet area, designed to measure susceptibility of a jet to soft radiation. In a similar 
manner, a quantity measuring the intra-event fluctuations, denoted as a, is determined from the 
sorted list {ptj/Aj}. It is given by the value for which 15.86% of jets have smaller ptj/Aj. 

3 Understanding systematic effects: a toy model study 

It is not guaranteed a priori that the methods described above will give a sensible result for 
the characteristic momentum scale of UE. Therefore, to better understand how those methods 
perform, we have tested them against a simple toy model. 

The model involves two components: a soft one, which we identify with genuine UE, and 
a hard one, which comes from perturbative contamination. The main parameters of the soft 
component are: the average density of particles per unit area, v, and the average pt of a particle, 
p. Therefore, the true value of p = vp by definition. The number of particles is governed by the 
Poisson distribution and the pt of a single particle by the exponential distribution. The hard 
contamination is modeled as coming from the initial state radiation of soft and collinear primary 
emissions. We assume this radiation to be independent of rapidity and <fi for 1 < pt < Q, where 
Q is a hard scale of the process (e.g. half of the hardest jet's pt for the dijet event). 

One of the results of the toy model study are the biases for p and a extraction. ^ For 
instance, the hard radiation introduces the bias for p determined in the traditional approach, 
bp ~ c? s Q, which depends linearly on the hard scale of the process. The corresponding bias 
in the area/median approach has only log(log) dependence 5p ~ a (rib/A tot + const • log log Q), 
where rif, is the number of final state born particles and Atot is the total used area of the event. 

Another important outcome of the study with the toy model concerns fluctuations. The 
extracted values of p vary from one event to another, even if the same p distribution is used 
to generate all events. This is because, usually, one works only with limited part of an event 
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Figure 1: Left: Distribution of p e xt from toy model. Right: Rapidity dependence of (p) from series of MC model. 



(y and eft cuts) and a finite number of objects (particles, jets). The magnitude of these intrinsic 
event-to-event fluctuations is an important characteristic of a method. 

Fig. Q] (left) shows the histograms of p from the toy dijet events extracted with the traditional 
method in its TransMin and TransAv variants and with the area/median method. The latter 
was used after removing the two hardest jets from the list of jets. We see that using the 
area/median based method results in the peak which is better both in terms of position and 
the width. The real difference is seen, however, in the values of the standard deviation. In the 
traditional method, it comes out as big as the value of p itself, while staying moderately small for 
the area/median method. The reason for this is again related to the difference in the sensitivity 
to the hard radiation between average and median (linear vs log(log) as discussed above for p). 



4 Approaching real life: Monte Carlo study 

The weak sensitivity of the area/median method to the hard radiation makes it advantageous 
for event-by-event studies and for measuring fluctuations. Therefore, we have used this method 
to examine the underlying event from a series of Monte Carlo generators and tunes. ^ We have 
carried out the study for dijet events at y/s = 10 TeV. The leading and next to leading jets, 
found with the anti-fct algorithm^! with R = 0.6, were required to lie in the rapidity window 
\y\ < 4 and to have pt greater than 100 GeV and 80 GeV, respectively. For the study of the UE 
we used the C/A algorithm with R = 0.6. 

As a first step, we verified that the more realistic UE from MC models shows a number of 
characteristic features found in our study of the toy UE. Then, we have examined a series of 
observables of potential interest for measuring at the LHC. 

The main quantity, (p), is show in Fig. Q] (right) as a function of y. We note that the 
rapidity dependence is quite strong and the exact level of the UE extrapolated for the LHC 
depends on the MC model/tune. UE may fluctuate both from point to point within a single 
event and from one event to another. The first kind of fluctuation is measured by (<t) and is 
shown in Fig. [2] (left). We see that all the models predict large intra-event fluctuations. Another 
interesting thing to note is that in Herwig they are nearly 40% smaller than in Pythia. This 
difference is consistent with the correlation coefficients shown in Fig. [2] (right) as a function of y2 
for a single y\ bin. Altogether, the results from Figs. [Hand [2] illustrate the potential gain to be 
had from studying wider variety of observables. 
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Figure 2: Left: Average intra-event fluctuation (a) from series of MC models as function of rapidity. Right: 
Correlations of p(y2) with p(yi), as function of j/2 for y\ in the bin — 1 < y < 1. 



5 Conclusions 

We have carried out a twofold study devoted to the problem of measuring the underlying event. 
Using a simple toy model of UE, we have examined the methods of its determination. Subse- 
quently, by studying the events simulated with MC generators, we have identified a range of 
important characteristics of the UE. 

One conclusion from the toy model study is that for determinations of averaged quantities, 
like {p), both the traditional and the area/median measurement methods give comparably good 
results. In contrast, for event-by-event measurements and determinations of fluctuations of the 
soft component, the traditional approach is affected significantly more by the hard radiation. 

Therefore, we chose the area/median method to examine more realistic UE from the MC 
models. We found noticeable differences between predictions of different generators/tunes ex- 
trapolated to LHC energy. For this reason we advocate measuring a broader range of observables 
including rapidity dependence of p, intra- and inter-event fluctuations and correlations. 
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